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Abstract. Relativistic quantum plasma dispersion functions are defined and the 
longitudinal and transverse response functions for an electron (plus positron) gas are 
written in terms of them. The dispersion is separated into Landau-damping, pair- 
creation and dissipationless regimes. Explicit forms are given for the RQPDFs in the 
cases of a completely degenerate distribution and a nondegenerate thermal (Jiittner) 
distribution. Particular emphasis is placed on the relation between dissipation and 
dispersion, with the dissipation treated in terms of the imaginary parts of RQPDFs. 
Comparing the dissipation calculated in this way with the existing treatments leads to 
the identification of errors in the literature, which we correct. We also comment on a 
controversy as to whether the dispersion curves in a superdense plasma pass through 
the region where pair creation is allowed. 
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1. Introduction 

There is a relatively sparse but diverse body of literature on relativistic quantum 
plasma dispersion functions (RQPDFs). Tsytovich (1961) calculated the response of 
a relativistic quantum (RQ) electron gas and derived expressions for the real and 
imaginary parts of the longitudinal and transverse response functions for arbitrary 
isotropic distributions, and for a Jiittner distribution. Jancovici (1962) derived 
expressions for the real and imaginary parts of the longitudinal and transverse 
response functions for a completely degenerate electrons gas, providing a relativistic 
generalization of the well-known result of Lindhard (1954) for a nonrelativistic 
degenerate electron gas. These results were rederivcd and extended in various ways, 
by Hakim and Heyvaerts (1978, 1980) and Sivak (1985), using a Wigner function 
approach (Hakim 1978), by Delsante and Prankel (1980) and Kowalenko et al (1985), 
who concentrated on the longitudinal response, and by Hayes and Melrose (1984), 
who derived general results and Melrose and Hayes (1984), who extended Jancovici's 
(1962) results to the nearly degenerate hmit and Tsytovich's (1961) results to a mildly 
degenerate plasma. The response functions were discussed further by Itoh et al (1992) 
and Braaten and Scgcl (1993) in connection with neutrino losses from stellar interiors. 
This area is of ongoing interest (Ratkovic, Dutta and Prakash 2003; Dutta, Ratkovic 
and Prakash 2004; Kocrs and Wijers 2005; Jaikumar, Gale and Page 2005). One of the 
neutrino emission processes, the so-called plasma process, is dependent on dispersion and 
dissipation in the plasma, and the results of Braaten and Segel (1993) are used. Although 
these authors included RQ effects in their formal development, they approximated the 
general result, effectively by neglecting RQ effects in the resonant denominator, and 
it is their approximate form that is used. We argue that this leads to a misleading 
conclusion concerning one-photon pair creation (PC) in a superdense plasma, which 
we define here as plasmas with cutoff frequency (called the plasma frequency by some 
authors) exceeding the threshold, 2m, for PC. (We use natural units with h = c = 1.) 

Our main purpose in this paper is to discuss the properties of RQPDFs in an 
isotropic, thermal RQ electron gas emphasizing the relation between dispersion and 
dissipation and the role of PC. One general feature of existing treatments is that the 
dispersion and dissipation are treated in different ways, with the dispersion described 
by appropriate RQPDFs, and with the dissipation calculated directly from the resonant 
part of the response function. In principle, the RQPDFs have imaginary parts that 
describe the dissipation, with the relation between the real and imaginary parts 
determined by the causal condition (e.g., the Landau criterion). Calculation of the 
dissipative part in this way provides a useful consistency check on the expression for the 
real part of the response function. Our consistency check fails in two published cases, 
and we identify and correct the relevant errors. 

There are relatively few applications where the combination of intrinsically 
relativistic and quantum effects is important in a plasma. The extreme conditions 
required apply in, for example, the early Universe, quark-gluon plasma and the interiors 
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of compact stars. In section 6 we discuss possible applications, emphasizing a specific 
point relevant to the plasma process for neutrino emission: whether or not one-photon 
PC is possible is a superdense plasma. Earlier authors (Tsytovich 1961, Hayes and 
Melrose 1984, Kowalenko et al 1985) assumed that the dispersion curve does pass 
through the PC region, so that PC is allowed, and Itoh et al (1992) and Braaten 
and Segel (1993) gave arguments against this. Here we show that PC is possible in 
a superdense plasma and we determine the conditions under which it can occur. 

In section 2 we present general formulae for the response functions and define 
general forms of RQPDFs for isotropic distributions. In section 3 we discuss dissipation, 
giving particular emphasis to the boundaries of the regions when Landau damping (LD) 
and PC are allowed for a given particle. In sections 4 and 5 we discuss the responses 
of a completely degenerate distribution and a nondegenerate thermal distribution, 
respectively. We discuss applications in section 6 and summarize our conclusions in 
section 7. 

2. RQPDFs for isotropic distributions 

An isotropic distribution is isotropic in one inertial frame, which is the rest frame of the 
medium. The linear response tensor for an isotropic medium may be described in terms 
of the longitudinal and transverse response functions. In this section, we start with a 
general form for the response tensor, then write down explicit forms for the longitudinal 
and transverse response functions, and identify relevant RQPDFs. 

2.1. General form, for the response tensor 

The general expression for the response tensor has been written down in a variety 
of different forms. We start with a covariant form that is derived by analogy with 
the (unregularized) vacuum polarization tensor (e.g., Berestetskii et al 1971), V^^{k), 
which relates the Fourier transform in space and time of the linear induced 4-current, 
J'^{k), to the 4-potential, A'^{k), where k denotes the wave 4- vector, [a;,k], constructed 
from the frequency, uj, and the wave 3-vector, k. The response 4-tensor satisfies the 
charge-continuity and gauge-invariance relations, knV'^^{k) — 0, kj/P'^^{k) — 0. This 
form is 



where n''(p) is the occupation number, with = ±1 labeling electron and positron 
states, and with 




W - C) + C^^(P) - C^n^'(pO F^'^icp, CP') 



(1) 



u-Ce + C'e' CeC'e' 



F'^'iP, P') = PI'P"' + P'l'P" + gi"'{m'^ - PP'). 



(2) 
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The energies are e = e(p) = (m^ + |pp)^/^, e' = e{p'), and p denotes the 4-momentum 
with components [e,p]. The (unregularized) vacuum polarization tensor itself follows 
from (0) by neglecting the contribution of the particles (^^^(p) — > 0). 

2.2. Longitudinal and transverse response functions 

Isotropy implies that the response tensor is of the form 

V^'-'ik) = V^{k) L^'-'ik, u) + V^{k) T'"'{k, u), (3) 

where L'^'^{k, u) and T^'^{k, u) are longitudinal and transverse projection operators that 
depend on the 4-velocity of the rest frame of the plasma. After projecting to 
identify V^{k) and V^{k), one may perform the p' integral over the (5-function, choose 
the rest frame and rewrite the p-integral in terms of integrals over e, e': 

Sp^2nJ d|p| |p|2y ^dcos^ = — / dee I de'e', 



m J£ 



where the limiting values are 

The response functions are unchanged by interchanging electrons and positrons. One 
may evaluate the response functions for electrons, and then replace the occupation 
number by 

h{e) = n^{e)+n-{e), (5) 

to include the contribution of the positrons. 

Explicit forms for the response functions (neglecting the vacuum contribution) are 
(Hayes and Melrose 1984) 

^ ^ m|k|2 ^ 87r2|k|3 
V\k) 



e2npo(^^ + 


k 


2) e\uj^- 
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^)m 


2m k 
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3 



[uo^ - |kn5(°)(A;) - AmuoS'^^\k) + Am^S^^\k) 

44 + a;2 + 2|k|2)^(0)(fc) 
AmujS^^\k)+Am'^S^'^\k)\, 



(6) 



(7) 



with npo the proper number density, with 



|k| /cu2- |k|2-4m2y/' 

"'^ - T 1, - |kP ) ' 

and where three RQPDFs are introduced, S^"'\k), with n = 0,1,2. These RQPDFs 
involve integrals over the occupation number and two logarithmic functions: 

/ N r de / e \ "- f In Ai for = 0, 2, 
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The response functions V^'^ are related to the dielectric response functions e^'^ of 
Jancovici (1962) and Kowalenko et al (1985) as follows: 

Ime^'^ = ^ImP^'^, Ree^'^ = 1 + ^ReP^'^, (11) 

resulting in dispersion relations in the rest frame of the plasma of the form 

+ /ioReP^(fc) = 0, - |kp + /ioReP^(fc) = 0. (12) 

3. Dissipation 

As already noted, in the existing literature it has been conventional to treat dissipation 
separately from dispersion, rather than treating it in terms of the imaginary parts of 
RQPDFs. Dissipation is possible only when e^, as given by (jHl), is real: LD applies for 
uo'^ — |kp < and PC for uo'^ — |kp > Am?. Before discussing the conventional procedure, 
it is useful to identify the energy, momentum and speed of a resonant particle at the 
boundary of the allowed regions for LD or PC. 

3.1. Limiting values of the resonance condition 

In dispersion theory it is conventional to refer to a zero of the denominator (as a function 
of UJ for fixed k) as a resonance, and to the algebraic condition for such a zero as a 
resonance condition. The Landau prescription specifies how one is to integrate around 
the associated pole in the integrand in accord with the causal condition, such that each 
pole contributes an imaginary part equal to its semi-residue. Dissipation is described by 
the contributions from these semi-residues, and different poles are interpreted in terms 
of different dissipation processes. There are two dissipation processes in a collisionless, 
unmagnetized plasma, LD and PC. 

For given cj, |k|, the limiting values of the resonance condition determine discrete 
values of £, |p|, |v|. These values are found by setting u) — C,e + ('e' = 0, squaring twice 
to remove the square roots, and setting (p ■ k)^ = |pp|kp. Writing |p| = msinhx, 
£ = mcoshx, |v| = tanhx, t = tanh(|x), the limiting values must satsify 

u;^ - IkP Ikl 

(l + a)¥ + 2(l-a^-262)t^ + (l-a)2 = 0, a= —, b = —. (13) 

The solutions for t^ are 

e = 4, t± = — — '—. 14 

1 + a 

There is considerable freedom in choosing the four solutions. An obvious choice is 
t = ±t_|_, t = ±t_. Noting that t and —1/t correspond to the same values of e, |p|, |v|. 
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one may also choose from t = ±l/t+, t = ±l/t_. With t = t±, the solutions for the 
energy, momentum and speed are 

e± _ a± b{a^ + b^- 1)^/^ p± _ ab ± {a^ + - 1)^/^ 
m 1 — 6^ ' m 1 — 6^ ' 

p± b±a{a' + b^-lY/^ . 
v± = — = i^—r^ • 15 

Another form of the boundary solutions was introduced by Tsytovich (1961); these are 
related to the i-solutions by 

e± = i^u±ek, p± = ^\k\±—ek, (16) 

with Ek given by (jH)). The solutions (fT!)|) are the natural solutions for PC, and they 
correspond to the choice t = t±. For LD, the natural solutions are e = ± 
IpI = c<j£:fe/|k| ± ||k|, and these correspond to t = t_ and t = l/t+, respectively. At a 
boundary of the LD or PC region, e must correspond to either e±, with e'maxmin then 
corresponding to e^, but some care is required in making the specific identifications. 



3.2. Dissipation due to LD and PC 

The resonant terms in (Q) correspond to replacing the denominator by —i7iS{uj—(e+('e'), 
where the Landau prescription is used. We discuss LD and PC separately. 

The terms with C = C = describe LD, and these correspond to resonances 
at a; = ±{e — e'), respectively. On repeating the derivation of the longitudinal and 
transverse parts, the double integral over e, e' is reduced to a single integral by this 
(5-function, and the limit of integration can be expressed in terms of Sk- In this way 
Tsytovich (1961) derived imaginary parts that correspond to 

p^i.p' roc 

ImVUk) = ^ de" {As'" - \k\')[n{e" - ^uj) - n{e" + Ic.)], (17) 

IraVlM = ^'g|n^|3 4 d." {2e"' - 2el + ^[^6" - ^u) - n{e" + i^)], (18) 

in the notation used in this paper. 

The resonant part of describes PC for CQ' = —1, but only ( = 1 = — C 
contributes for uj > 0. The vacuum contributes in this case, and Tsytovich (1961) 
retained both the contribution of the vacuum and of the electron gas, but his final 
expression for the vacuum contribution differs from the well-known result by a factor 
of two. We note that the response is invariant under the interchange of positrons and 
electrons, which allows one to rewrite 1 — n'^{e) — n~{e') in for ( = 1 = as 
1 — ^n{e) — |n(e'). Then, repeating the derivation, we find 

ImV^cik) = £1 de" {As"' - ^[1 - + e") - \n{\u - e'% (19) 

^^ndk) = ~ '^''^ r d5"(25"2-24 + |k|2)[l-in(iu; + .")-Hiu;-5")],(20) 
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for the imaginary parts due to PC of the longitudinal and transverse responses, 
respectively. The unit term inside the square brackets differs by the relevant factor 
of two from Tsytovich's expression. 

3. 3. Dissipation due to the vacuum polarization tensor 

The vacuum polarization tensor is of the form Pq''(A;) = Vo{k){g^'^ — k'^k'^/k'^), 
where Vo{k) is an invariant. The longitudinal and transverse parts are V^lk) = 
(jj'^Voik) / {u}"^ — |kp) and Vji^k) = Vo{k), respectively. The real part is negligible for 
most purposes involving wave dispersion, but the imaginary part cannot be neglected 
when considering dissipation due to PC. An important point is that dissipation due to 
PC occurs in the vacuum, and the presence of an electron gas tends to suppress it due 
to the Pauli exclusion principle. The imaginary part of the vacuum polarization tensor 
is well-known, e.g., Berestetskii et al (1971), and corresponds to 



for io'^ — |kp > 4m^, with ImPo(^) = for lo'^ — |kp < 4m^. A derivation of (j2II) using 
the approach adopted here leads to (fT^ and (pUj) with only the unit terms retained. 
The integrals are then elementary and (fTT^ and reproduces the longitudinal and 
transverse parts of (PT|) . respectivley. 

3.4- Imaginary parts of logarithmic PDFs 

Logarithmic functions appear naturally in Q, and more specifically as RQPDFs for 
a completely degenerate electron gas, as discussed below. It is desirable to have a 
prescription that allows one to write down the imaginary part of a logarithmic PDF 
directly. Superficially, this seems trivial: when x passes from positive to negative, \n.x 
may be replaced by In ± ivr. However, determining the relevant sign is not trivial. 

The imaginary part of any PDF may be determined by imposing the Landau 
prescription. For a logarithmic PDF, this leads to the generic prescription 



Hence, to impose the causal condition, one writes a logarithmic function as a sum (or 
difference) of terms of the form (j22j) and gives each term an imaginary part, ivr, when 
the frequency is in the range uJunn < uJ < cUmax- 

4. Completely degenerate Fermi gas 

The response of a completely degenerate Fermi gas was calculated by Jancovici (1962), 
cf. also Hayes and Melrose (1984), Sivak (1985), Kowalenko et al (1985). Jancovici's 
expression for the transverse part of the response tensor contains a spurious factor 



l^Vo{k) = ^[m' + \{u'-\\.n^ 

OTT K 



(21) 




max; 



(22) 
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uj^jiijp' — |kp), which leads to a nonphysical resonance in the dispersion relation for 
transverse waves. It also implies incorrectly that the contribution of the electron gas 
to dissipation due to LD and PC has the same sign for the transverse response. Here 
we start with forms that are valid in the DL region, where the response functions are 
necessarily real, and then discuss the extensions into the LD and PC regions. 



J^..^. Thermal distributions 

Before considering the completely degenerate limit, it is appropriate to comment on the 
general case of a thermal distribution of electrons, which is the Fermi-Dirac distribution, 

""^^^ " exp[{e-fi)/T] + l + exp[{e + fi)/T] + V ^^^^ 
where the chemical potential, fi, includes the rest energy, m, and where the temperature 
is T, with Boltzmann's constant set to unity. (The form ()23|) applies for each spin state, 
and a factor of two arises when one sums over the two spin states for unpolarized 
electrons and positrons.) No general results are known for the RQPDFs for this case. 
The completely degenerate limit corresponds to T —>■ 0, fi ^ ep, when one has 

{1 for e < ep, 
(24) 
for e > ep, 

where ep = {m? +p|)^/^ is the Fermi energy, with the Fermi momentum determined by 
the electron number density, rie =pp/3n'^. 



4-2. DL region 



In the DL region Jancovici's response functions, in the present notation (with the 
spurious factor omitted), are 



e^u^ (Seppp 2\k\'^ f ep + pp\ , ep[Ael + 3{uj^ - 



47r2|k|2 I 3 
c^[3|k|2 



+ 



6 k 



In A 



IF 



+ 



UJ 



47r2 



12|k| 

4cu2 + 2 |kp 
3|k 

+ep 



■ eppp + 



InAaF + 



2m2 + cu^ 




|k 


2 


3(cu2_ 


k| 





\k\ek -r-r InAsF 
\uj\ 



In 



Sp+pp 



m 



3|k| 



+ 



Am? 


k 


2 + 


k 


4 


4|k 


3 



In A 



IF 



c<;[12m2|k| 



|k|2)(124 + a;2 + 3|k|2)] 



2m2 + cj" 



3|k| 



24|k|3 

— £fc -; — r In A 



In A 



2F 



\L0\ 



3F 



with AjF given by setting |p| = pp, e = ep m. the expressions (fTUj) for Ai, A2 and 
As 





k 




k| + 2|p|£fc)2 - 4m^tj2 


k 


2 


(Cj2 - 


k 


me 


k 


-2|p|£fc)2 -4m4cj2 


k 


2 



(25) 



(26) 



(27) 
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Note that the sign t^/jco'l is needed in the terms involving A3F to ensure that the real 
parts of the response functions are even functions of u. 

The DL regime corresponds to |kp < u"^ < Am? + |kp. In this regime, Aip and 
A2F are positive (they are always real) so their logarithms are real. Due to Sk being 
imaginary, A3F is the ratio of a complex number and its complex conjugate, and hence 
its logarithm is imaginary, so that e^lnAsF is real. One may write 

, arctan y 1 — arctan y_ for y_ > 0, 

arctan x+ — arctan X- — tt for X- < 0? 



X± = i,,iro,. oii.i , „ ro iTJ^- (29) 



2£f(^^ - 


|k| 






|k|[2m2|k| ±pf{l 


^2- 


|kP)] 



Although earlier authors, e.g., Kowalenko et al (1985), Sivak (1985), noted that In A3F is 
replaced by an arctangent, there are many choices as to how this arctangent is written. 
With the choice made in ()28p the arctangents x± remain between and 7r/2 throughout 
the DL range. The choice avoids complications with other choices in numerical 
calculations. 

4.3. Imaginary parts o/lnAjF 

In order to use the prescription (j221), the logarithmic functions must be written in an 
appropriate form. Relevant forms for Ai, A2 follow by writing ()10j) in terms of the 
limiting values (j^. This gives 

^ ^ (^ - g + g^ax) (^ + g - gmax) (^ - g - ^max) (^ + g + ^max) /ggN 

^ {uj-e + e'^^^) {u; + e- £'„,;„) {uj - e - e[^^^) {u; + e + e'^^^) ' 

^ (^ - g + g^ax) (^ - g - gmax) (^ + g " gmin) (^ + g + gmin) j^g^^. 

Although A3 cannot be rewritten in terms of the factors that appear in ()30|) and (jHTj) . 
it can be written in a form similar to (jHT|) . with uj replaced by 2ek'- 

j^^ ^ (2gfc - g + gmax)(2gfc - g - gmax)(2gfc + g - gmin)(2gfc + g + gmin) ^^2) 
(2efc - g + gmin) (2gA: - g - gmin) (2gfc + g - gmax) (2gA: + g + gmax) 

The logarithms of Ai, A2 may be written as a sum of terms of the form by identifying 
(a;-a;max)/(^-Wmin) with {Lj±e s'^^^) / {LJ±e ~ e'^^J or with {uj±e+e'^iJ/{LJ±e+e'^^^). 
Although In A3 cannot be rewritten in terms of these factors, in the neighborhood of 
the zeros of any of the factors in (jH^ . the vanishing factor does become of this form. 
As only the sign of the imaginary part on crossing the zero is required, this suffices to 
determine the sign. For example, consider the factor {26k — s — gmax)- ^^is factor is 
zero at the boundary of the LD region with e = Ek — gmax = ^fc + and in the 
neighborhood of this boundary the factor may be approximated by {u — e — g'max) ' 
treated in the same manner as the corresponding factor in Ai or A2. 

The boundaries of the allowed regions for LD and PC are illustrated in figure Q 
For |k| > 2pF, one has ep < gprnini then the upper and lower frequency boundaries 
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Figure 1. Regions of w-|k| space {u > 0) are separated by curves corresponding 
to tlic boundaries of the regions wliere LD is allowed, uj < |k|, and PC is allowed, 
uj > (4m2 + |kp)i/2. These are further separated into regions (a)-(h) defined for a 
completely degenerate electron gas with p-p/ni = 1.5. For the completely degenerate 
gas, LD is allowed only in regions (b) and (c), and there is no dissipation in (a) 
and (d), the electron gas completely suppresses PC in (e) and partly suppresses 
PC in (f); PC has its vacuum value in (g) and (h). From lower right to upper 
left the curves are: uj = (el — 2pF|k| + |kp)^/^ — ep (dashed, solid diamonds), 
UJ eF-(e|-2pF|k| + |k|2)i/2 (dashed, solid squares), to = (e^ + 2pF|k| + |k|2) V2 _£p 
(dot-dashed, open squares), uj = |k| (sohd, solid circles), us = {Atit? + |k|^)^/^ (dotted, 
open circles), w = (e| — 2pF|k| + |k|^)^/^ + ep (double-dot-dashed, open triangles), 
UJ — {ep + 2pF|k| + |k|^)^/^ -I- ef (dashed, solid triangles). 

Table 1. Imaginary parts of RQPDFs for a completely degenerate electron gas. 

upper boundaries In Ai In A2 In A3 

(c) LD UJ < Ef - epmin^ |k| < 2pF 

(b) LD leF-^minl <^<4max-^F IVT 

(e) PC (4m2 + |k|2)i/2 < < £f + ^Fmin, |k| < 2pF 

(f) PC ep + e'p^i^ <u <eF + e'p^^^ m 

are u = ^Fmax.min ~ LD, and UJ = ^Fmax.min + for PC. In this case only 

ln[(ti; ± £:f — £^Fmax)/(^ ± £:f — f^Fmin)] Contribute to LD and PC, respectively, with these 
factors giving an imaginary part of ivr in the region where the argument of the logarithm 
is negative, and zero otherwise. For |k| < 2pF one has ^f > ^Fmin- ^^i^ case, the zero 
of — ^F + £^Fmin occurs withiu the LD region, separating regions (b) and (c) in figure UJ 
and the zero of uj — ep — ^Fmin occurs within the PC region, separating regions (e) and 
(f) in figure n It is then straightforward to determine the signs of the imaginary parts 
in the various regions, and these are listed in table 1. 

The imaginary parts oiV^{k) and V'^{k) may be written down by inspection using 
()25|1 and ()2(j|l . respectively, and noting the imaginary parts in tabled These imaginary 
parts may also be derived from (fTTj) and (fTHj) by setting the occupation number equal 
to unity for e < e-p and zero for e > e-p, and performing the integrals, which are then 
elementary. The imaginary part was written down for the longitudinal response by 
Jancovici (1962) and Kowalenko et al (1985) and these results are reproduced using 
^ and tabled 

4.4- Alternative forms for Ki 

For completeness we note that the logarithmic functions may be written in terms of the 
± solutions (HU) and (fT3)): 

^ ^ (t + t+)(t + t-)(t-l/t+)(t-l/t_) _ (|p|+P+)(|p|+P-) 
1 (t - t+){t - t^){t l/t+){t + l/t_) (|p| - p+)(|p| - p_) ' 



-1271 

—ivr ivr 

-2i7r 

ivr —ivr 
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(|v| 


+ '^+)(|v| 


+ V-) 


(|v| 


-^+)(|v| 




(|v| 


+ ^+)(|v| 




(|v| 


-^+)(|v| 


+ v^) 



^ _ (t + t+)(t + t_)(t + l/t+)(t + l/t,) 

The boundaries in figure ^ are identified as follows: the upper boundary to the PC 
region corresponds to tp = —t- (pp = —P-, Vf = —V-), the lower boundary of region 
(f) corresponds to tp = t_ {pp = p_, vp = f_) for |k| > 2pp and to tp = t+ {pp = p+, 
Vp = v^) for |k| < 2pp; the upper boundary of region (b) corresponds to tp = 1/^+ 
{Pp = —p+, Vp = 17+ ), the lower boundary corresponds to tp = — 1/^+ (pf = P+, 
Vp = — f+) for |k| > 2pp and to tp = t_ {pp = Vp = vJ) for |k| < 2pp. Although 
useful for some other purposes, these forms are not convenient for determining the 
imaginary parts because the frequency-dependence is implicit rather than explicit, and 
the prescription (j22|l cannot be used directly. 



5. Nondegenerate thermal distribution 

The nondegenerate limit of the Fermi-Dirac distribution ()2H|1 applies when /z/T is large 
and negative, and then it becomes the Jiittner distribution, 

,X,, = ^e-'/-, A = 2 cosh(,/r) = (34) 

where K2 is a modified Bessel function, and he is the number density of electrons plus 
positrons. (The normalization coefficient, A, is evaluated by setting the integral of h^e) 
over 2d^p/(27r)^ equal to fig, where the factor 2 arises from the sum over the two spin 
states.) 

According to Melrose and Hayes (1984), in this case the three plasma dispersion 
functions S^^'^k) can be evaluated in terms of the RPDF introduced by Godfrey et al 
(1975), which they wrote in the form 

/■I df 

T{z,p)= exp(-p7), (35) 

J-i V — z 

where z = C(j/|k| is the phase speed, 7 = (1 — v"^)^^^"^, and p = m/T is an 
inverse temperature in units of the rest energy of the electron (p = 1 corresponds 
to T = 0.5 X 10^° K). Dissipation is described by the imaginary part of this RPDF: 

f ne-f^° for \z\ < 1, 
lmT{z,p) = { (36) 
[ for \z\ > 1, 

with 7o = (1 - ;z2)-i/2_ The RQPDFs become 

P ± l±v± VP / 



5«(fc) = ^E 
p ± 



P v± \ p 
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S^'^ (A:) = - E ^1 (4 + 7l) f ^ T'K, p) + 2K,ip)) 

P ± 'y±v± I VP / V P / 

- 2-fivlK,{p) - ^ 7l^±[r(t;±, p) + 2vMp)]], (37) 

with T'(t;±,p) = dT{v±,p)/dv±, 7± = (1 - and with A given by The 

sign a± = e±/\e±\ is needed in the LD region to take account of the fact that e± can 
be negative while 7± = \e±\/m are positive by definition. The sign a± is replaced by 
unity in the DL and PC regions. The RQPDFs ()H7|) characterize the response of a 
nondegenerate thermal electron gas when RQ effects are included. 

The interpretation of v± and 7± is different in the LD and PC regimes. In the LD 
region one has 

^± = rM Z — 7^' l± = {ek±2^)/m, (38) 

|K| Ek ± iO / Z 

with Sk defined by (jHI)- In the nonquantum limit, one has Sk ^ Ci^/2, |kp/2c<j, implying 
v± — i> L(j/|k| = z. Then fl37|) reproduces the known non-quantum limit (Melrose and 
Hayes 1984): 

S'^'\k) = '^J,zT{z,p) + 2K,{p)l 
2A 

5(1) (fc) = — [ - zT{z, p) + (1 - z') T'{z, p) + 2pK,{p% 
P ^ 

S^'\k) = ^[^^^zT{z,p) + 27o^(l + z')Mp) + K,ip)]. (39) 

In comparing (jHTj) with (jH^jl . it is apparent that the phase speed, z = a;/|k|, in the 
nonquantum case is replaced by two functions v± that include the effect of the quantum 
recoil, which has opposite signs for emission and absorption. Thus, in the LD region, v± 
are interpreted as resonant phase speeds for stimulated emission and true absorption, 
which differ due to the quantum recoil, and m7± are interpreted as the energies of the 
electron before and after emission of a wave quantum, respectively. In the PC regime, 
m7± = Ek ^ ^uj are interpreted as the energies of the created electron and positron. 

In the DL regime the RQPDFs must be real. In this case, the v± are complex 
conjugates of each other. With a± = 1 and T{v*,p) = T*{v,p), the sum over ± in (j?fj) 
leads to a real expression, as required. 

Comparison of the imaginary parts for the response functions obtained from the 
imaginary parts of the RQPDFs with those obtained by imposing the causal condition 
directly provides a check on both results. Here the imaginary parts that are to be 
compared are those obtained by inserting the imaginary parts of the S^'^\k) into the 
expressions (jH)) and ((Tj) for the longitudinal and transverse response, and those obtained 
by evaluating the integrals in (fTTj) . (fTHj) and (fT^ . (pH) for the nondegenerate distribution 
The results agree. In particular, Tystovich (1961) wrote down explict expressions 
for dissipation due to LD and PC in a nondegenerate electron gas, which our derivation 
reproduces: the spurious factor of two is present only in Tsytovich's calculation of the 
vacuum contribution. 
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6. Applications 

RQ effects become important for dissipation and dispersion in plasmas only under 
extreme conditions, such as the early Universe, quark-gluon plasmas and the interiors 
of compact stars. The degeneracy condition (temperature less than chemical potential) 
is relevant to only the last of these, and we concentrate on this case, emphasizing the 
role of PC. First we comment on LD. 

LD in the non-quantum limit is possible only for subluminal waves, u; < |k|, and 
this is also the case when RQ effects are included. This precludes LD of transverse waves, 
which are superluminal. Dissipation and dispersion associated with LD may be treated 
nonrelativistic provided not only that the particles are nonrelativistic, but also that the 
waves are subluminal. The RQ recoil term changes the classical resonance condition to 
u; — k-v±(ci;^ — |kp)/2m7, for emission and absorption, whereas in nonrelativistic theory, 
the recoil term is =F|kp/2m. The difference is unimportant for nonrelativistic particles 
and subluminal waves, but is important for superluminal waves and for waves with near 
vacuum dispersion, u ^ |k|. An implication is that the widely used response functions 
of Lindhard (1954), which were derived using nonrelativistic quantum mechanics, may 
lead to unreliable results for waves with > |k|. We are currently investigating this 
point. 

The contribution of the plasma process to neutrino emission from the cores of 
compact stars depends on the dispersive properties of the degenerate gas, and Braaten 
(1991) pointed out that earlier authors, following Baudet, Petrosian and Salpeter (1967), 
had used inaccurate forms for the dispersion relations. A controversial point was raised 
by Braaten (1991), who criticized the claim by Baudet, Petrosian and Salpeter (1967), 
and subsequent authors, that PC needs to be taken into account in sufficiently hot and 
dense plasmas: if PC is allowed then photons decay into pairs much faster than they 
would decay into neutrinos. This point was discussed further by Itoh et al (1992) and 
Braaten and Segel (1993), who also concluded that PC is forbidden in a completely 
degenerate electron gas. In more recent discussions (Ratkovic, Dutta and Prakash 2003; 
Dutta, Ratkovic and Prakash 2004; Koers and Wijers 2005; Jaikumar, Gale and Page 
2005), the approximations of Braaten and Segel (1993) to the dispersion functions have 
been used. Our results show that PC is allowed in a completely degenerate electron 
gas, and the reason that our results differ from those of these earlier authors can be 
understood as follows. 

Itoh et al (1992) considered wave quanta at the cutoff frequency, cjc, and argued 
that although one can have uJc > 2m in a superdense plasma, the actual threshold for 
PC is 2eF, and one cannot have Uc > ^ep. The higher threshold is because all electron 
states below the Fermi energy are occupied, and for |k| = the electron and positron 
energies are equal. This argument is consistent with our results, but it applies only 
at |k| = 0. PC is forbidden in region (e) in figure ^ and this region shrinks as |k| 
increases. A dispersion curve for transverse waves that starts at uJc > 2m for |k| = 
is necessarily in region (e) for sufficiently small |k|, but then necessarily enters region 
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(f), where PC is allowed, before approaching the light line asymptotically. Braaten and 
Segel (1993) made approximations to the wave dispersion by neglecting the quantum 
recoil: specifically, if one combines the denominators in (0), the common factor may be 
written as (fc-u)^ — (fc^/2m)^, and Braaten and Segel (1993) argued that for practical 
purposes one can neglect the (/c^/2m)^ term. However, near the cutoff |k| — > 0, one has 
{ku)^ — (fc^/2m)^ —> a;^7^(l — a;^/4m^7^), and their approximation requires Uc <^ 2m7, 
which is not satisfied for all 1 < 7 < ^p/^^ in a superdense plasma. This approximation 
effectively excludes dispersion due to PC, and it is inconsistent to use it to argue that 
PC cannot occur. Our results show that PC does occur over a limited range of |k| in 
a completely degenerate electron gas, but not in region (e) in figure HJ due to exact 
cancellation of the vacuum contribution to PC by the electron gas. For a partially 
degenerate electron gas, even in region (e) the cancellation is not exact, and PC occurs. 

The original argument of Braaten (1991) against PC was based on mass 
renormalization of the electron, arguing that this suppresses the cutoff frequency, 
keeping it below the PC threshold. We do not comment specifically on this argument 
here. Our conclusion is that the arguments by Itoh et al (1992) and Braaten and 
Segel (1993) that neglected mass renormalization do not negate the original argument 
of Baudet, Petrosian and Salpeter (1967) that PC needs to be taken into account when 
considering the plasma process for neutrino emission in a superdense plasma. 

7. Conclusions 

In this paper we discuss the properties of RQPDFs for an isotropic, unmagnetized 
plasma. The dispersion is related to the dissipation, which includes the familiar Landau 
damping (LD), modified by the quantum recoil, and one-photon pair creation (PC). 
It is necessary to treat the dissipation in the LD and PC regimes differently, and to 
interpret them differently. LD has the same interpretation as in a non-quantum plasma, 
except that the resonance at the phase speed, z = uj/\^\, is replaced by resonances at 
two speeds, v±, and corresponding energies, e±, given by (IHHj) and interpreted as the 
resonant values for induced emission and true absorption when the quantum recoil is 
included. Dissipation due to PC in the electron gas has the opposite sign to LD and 
a different interpretation: PC exists in the vacuum, due to the imaginary part of the 
vacuum polarization tensor (j2II), and the presence of an electron gas partly suppresses 
PC due to the Pauli exclusion principle. 

An objective in this paper is to relate dissipation and dispersion by deriving the 
dissipation from the imaginary parts of the RQPDFs, which requires that the imaginary 
parts be determined explicitly. In particular, the logarithmic RQPDFs that appear 
for a completely degenerate electron gas acquire an imaginary part of iivr when their 
arguments becomes negative, and a prescription is needed to determine the sign of 
this imaginary part uniquely. We start from the DL region, where the imaginary part is 
necessarily zero, and analytically continue into the regions where LD and PC are allowed. 
We show that the Landau prescription leads to a relatively simple prescription for 
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identifying the sign of the imaginary part acquired when the argument of the logarithm 
changes sign. We compare our results with existing expressions for the imaginary 
parts derived in other ways, and find agreement provided that some minor errors are 
corrected. Specific errors identified are a spurious multiplicative factor in Jancovici's 
(1962) transverse response function and a factor of two in the expression derived by 
Tsytovich (1961) for the vacuum contribution to dissipation due to PC. 

In the absence of any plasma, dissipation due to PC is determined by the imaginary 
part of the vacuum polarization tensor for u; > {rn^ + |kp)^/^, and is zero otherwise. 
The presence of an electron gas tends to suppress PC, and the presence of a completely 
degenerate electron gas can completely suppress PC. Complete suppression at a given 
uj, |k| occurs if all potential states for the created electron are below the Fermi level. 
Although it was argued by Itoh et al (1992) and Braaten and Segel (1993) that PC cannot 
occur in a superdcnsc plasma, where the cutoff frequency exceeds the PC threshold 2m, 
we show this is not the case for at least a range of |k| ^ 0. The earlier arguments of 
Baudet, Petrosian and Salpeter (1967) on the imphcations of PC remain vahd and need 
to be taken into account in detailed analyses (e.g., Jaikumar, Gale and Page 2005). 
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